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ABSTRACT 

We investigate the role of the environment on the colour and stellar population gradi- 
ents in a local sample of ~ 3500 central and ^1150 satellite SDSS early- type galaxies 
(ETGs). The environment is parameterized in terms of the number of satellite galaxies, 
Ng a i in each group. For central galaxies, we find that both optical colour and mass- 
to-light (M/L) ratio gradients are shallower in central galaxies residing in denser 
environments (higher N ga i). This trend is driven by metallicity gradients, while age 
gradients appear to be less dependent on the environment and to have a larger scatter. 
On the other hand, satellites do not show any differences in terms of the environment. 
The same results are found if galaxies are classified by central age, and both central 
and satellite galaxies have shallower gradients if they are older and steeper gradients 
if younger, satellites being independent of ages. In central galaxies, we show that the 
observed trends can be explained with the occurrence of dry mergings, which are more 
numerous in denser environments and producing shallower colour gradients because 
of more uniform metallicity distributions due to the mixing of stellar populations, 
while no final clues about merging occurrence can be obtained for satellites. Finally 
we discuss all systematics on stellar population fitting and their impact on the final 
results. 
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1 INTRODUCTION 

Colour and stellar population gradients in g alaxies are pro- 
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trends with the mass have been recently investigated on 
samples of local ETGs, suggesting that different physical 
phenomena could concur to shape the gradients at low and 
high masses (e.g. T+ IO). From on e side , gas infa l l and 
supernovae feedback ([L arson 1974, 1975; Kawata 200 ll; 
Kawata fc Gibson] 20031; iKobavashll |2004| ; iPipino et al.l 



201(1 ITortora et aUhoilal ) seem to be the main phenomena 
driving the evolution of low mass ETG s , while merg- 
ing and AGN feedback l|Kobavashil |2004| ; iHopkins et all 
120091 ; ISiiacki et al.1 l2007h would work at larger masses 
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l|Dekel fc Birnboim| [2006). However, environment also plays 
a crucial role in galaxy evolution, since many physical 
phenomena like tidal interactions, strangulations and 
harassment would affect the star fo rmation in low mass 
galaxies (e.g. IWeinmann et al~ll2009h and ca use inner pop- 
ulation gradients at different mass scales (e.g. ITortora et all 
l2011al h In more massive systems, though, the major player 
in driving the ste l lar population mixing is galaxy mergin g 
lIDavis et al.l 1 19851 ; ISpringel et aill2005l ; iRomeo et ai1l2008f ). 
iKobavashil (|2004l ) have shown that overall merging events 
tend to flatten the metallicity gradients with time. Look- 
ing in more details, merging events can have a complex 
taxonomy: thus, minor or major merging, or even gas-rich 
or -poor merging are expected to produce different stellar 
population gradients. For instance, after the initial gas rich- 
merging events (generally occurring at high redshift), larger 
central metallicity and positive age gradients are observed 
|Kobavashil 12004 iMihos fc Hernquistl Il994h ; subsequent 
gas poor-merging may dilute the positive age gradient with 
time as well as make the metallic i ty gradients flattened out 
l|Whitelll98"ol . lHopkins et al.ll2009l . lDi Matteo et al.ll2009l ). 
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Although higher mass systems would be the ones which 
have experienced a larger fraction of merging events, it is 
also interesting to investigate if central galaxies in groups 
and clusters may differ, in their population gradient prop- 
erties, from satellite systems orbiting in the cluster poten- 
tial. The central galaxies in clusters (and groups) are the 
most luminous and massive (in terms of stellar mass) ob- 
jects in the Universe. They are found to have different lumi- 
nosity profiles wh en compared with typical cluster galaxies 
l|Schombertlll986l l and do not seem to be drawn from the 
same luminosity function as bright ellipticals (|Dresslerlll97sl . 
iBernstein fc Bhavsarl l200lh . These evidences suggest that 
the evolutive processes of central galaxies can be strongly 
different from the ones driving normal (satellite) systems. 
Moreover, they still hold imprints from the early evolution- 
ary stages since they reside in the very central regions of 
clusters and groups, where mass started to be accreted ear- 
lier after the Big Bang than other density environments. 
These regions have witnessed during the cosmic time a 
variety of galaxy i nteractions with th e environment, and 
merging events (e.g. iRomeo et al.| [2008) contributing to the 
mass accretion of larger and larger galaxy systems, a pro- 
cess which is more efficie nt in denser environments like 
group/cluster haloes (e.g. IWhilev et al. 1 120081 . IStott et al.l 
120081 IStott et al.ll2010l i 

From this perspective it is natural to expect that the 
galaxy in the center of very dense environments might be 
more sensitive to the effect of the large amount of merging 
events expected in the hierarchical growth which shall be 
recorded in the stellar population parameters (e.g. age and 
metallicity). 

Thus, we have consid ered a local sample of SDSS galax- 
ies (|Blanton et al.ll200 5hj). recently analyzed and discussed 
in T+10 and T+H, where we have discussed colour, mass- 
to-light ratio and stellar population gradients in terms of 
mass and compared with independent observations. These 
results have been framed within the predictions of hydro- 
dynamical and chemo-dynamical simulations of galaxy for- 
mation. In the present paper we will discuss the connection 
with the environment, selecting those galaxies classified as 
centrals and satellites in groups and clusters and investigat- 
ing if gradients change as a function of the environment. 

The paper is organized as follows. The data sample and 
the analysis have been presented in Sj2j the results are dis- 
cussed in 33 while 3H is devoted to the physical interpre- 
tation and conclusions. Systematics in stellar population fit 
have been discussed in iJX] 



2 DATA AND SPECTRAL ANALYSIS 

We start from a database consisting of 50 000 low redshift 
(0.0033 2 sC 0.05) galaxies in the NYU Value-Added 
Galax y Catalog extracted from SDSS DR4 (|Blanton et al.l 
12005b], hereafter B05jE recently analyzed in T+10 and 
T+llQ where stellar population synthesis models have been 
used to determine galaxy stellar mass, M*, colour and stellar 



population parameters/gradients. Following T+10, we have 
sorted out ETGs by keeping those systems with a Sersic 
index satisfying the condition 2.5 ^ n ^ 5.5 and with a con- 
centration index C > 2.6. The final ETG sample consists 
of 10 508 galaxies. We have cross-matched our datasample 
with the (z > 0.0 1) DR4 SDSS based group catalog from 
lYang et al l (|2007r ) to recover information about the envi- 
ronment the galaxies live. lYang et al l (|2007T ) have identified 
the groups and have separated the most massive (or most 
luminous) galaxies in each group, labeled as centrals, from 
the satellites. We will use both 1) the identification of satel- 
lite and centrals on the base of stellar mass selection, and 
2) the number of galaxies iVg a i in each group as an environ- 
ment indicator. We have retained galaxies - isolated (having 
Ngzi = 1), central and satellite - with a mass logM, > 10.5 
and left with a final sample of 3525 central (including 1941 
isolated systems) and 1141 satellite galaxies. As a further 
criterion, we have also ranked satellites in each group on the 
basis of their stellar masses. 

We notice that the fraction of central ETGs in each 
mass bin is lower than the fraction of satellite and field galax- 
ies at log M»/Mq < 10.7 — 10.9, while it increases at larger 
masses, with 62 centrals, 6 satellites and 1 field galaxy for 
log M*/M > 11.4. 

As discussed in T+10, we have used the structural pa- 
rameters given by B05 to derive the color profile (X — Y)(R) 
of each galaxy as the differences between the (logarithmic) 
surface brightness measurements in the two bands, X and 
Y. The stellar population p r opert ies are derived by the fit- 
ting of iBruzual fc Chariot I ^003, hereafter BC03) " single 
burst" synthetic stellar models to the observed colours . Age 
and metallicity are free to vary, and a lChabrier I (|200lh IMF 
is assumed. However, in order to check the effect of the ex- 
isting degeneracies between age an d metallicity, we will also 
assume the age gradient to zero (|Tortora et al.ll2011bl ) by 
fixing the age to 10 Gyr as a prioi0. We define the CG as 
the angular coefficient of the relation X — Y vs log i?/i? c ff , 

r/ v" V\ 

Vx_y = — — ttttt: — c> measured in mag/dex (omitted in 
Slog(R/R eS ) 

the following), where R a g is the r-band effective radius. By 
definition, a positive CG, Vx-y > 0, means that a galaxy 
is redder as R increases, while it is bluer outward for a neg- 
ative gradient. The fit of synthetic colours is performed on 
the colours at Ri = 7? o ff/10 and R2 = R c « and on the total 
integrated colours. Stellar parameter gradients are defined 
<51og(W0 

as Vw = — — , , _ , , where W = (ane,Z, T„) are the 
Slog(R/R Ba ) 

estimated age, metallicity and M/L. Because of the defini- 
tions adopted for Ri and R2, an easier and equivalent defi- 
nition can be used, in fact, we define V agc = log[age2/agei], 
V z = log[Z 2 /Zi] and V T » = log[T» 2 /T,i] where (age^, Z it 
T»i), with i = 1,2, are the estimated parameters at Ri and 
R2, respectively. 

In T+10 and T+ll we have discussed the impact of 
the wavelength coverage on the results. Here, we will add 
further Montercarlo simulations in order to verify that the 
optical bands alone introduce small spurious degeneracies 



1 The catalog is available at: 
http: / /sdss .physics .nyu. edu/vagc/lowz .html 

2 Details about sample selection, incompleteness and biases can 
be found in T+10. 



3 As we will also discuss later, V age = is a reasonable assump- 
tion for the older g alaxies, which are found to have V a g C ~ 
llTortora et aLlbOlOT) . 
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Figure 1. Gradients in terms of stellar mass for centrals (left) and satellites (right). In each panel we show g — r (top left), M/L (top 
right), age (bottom left) and Z (bottom right) gradients as a function of stellar mass. The medians and 25- 75th quantiles are shown. The 
dashed symbols are for the case with V agc = (i.e. agei = age2 = 10 Gyr). 



and/or correlations among parameters with respect to stel- 
lar models including near-IR or UV data. In particular, we 
will quantify any systematics and possible trends with stellar 
mass and N ga i in the Appendix I All and discuss the impact 
on our results in 



3 GRADIENTS AND ENVIRONMENT 

In T+10 and T+H we have discussed the colour, age, Z and 
M/L gradients as a function of stellar mass without distin- 
guishing central from satellite galaxies. Here, we will expand 
the analysis by investigating the effect of the environment 
on the most massive central and satellite galaxies, showing 
the significance of the correlations in Table [T] and the slopes 
in Tabled 

We start in Fig.[T]by showing the gradients V g _ r , Vt,, 
Vz and V age as a function of the stellar mass for central and 
satellite galaxies. The trend of central galaxies alone does 
not differ much from the one of the total sample studied in 
T+10 and T-fTl. They have colour and metallicity gradients 
which become shallower at larger masses. Vr.s are negative 
at low masses, turn out to be null around logAf*/M0 ~ 11 
and finally become slightly positive at very high masses. The 
positive trend of V age is statistically significant (as shown 
in Tables [T] and [5J) despite its large scatter. Vz show an 
increasing trend with M* (from -0.5 to -0.3 across the mass 
range), although the scatter at lower masses is quite large. 

If one forces V agc = (assuming agei = age^ = 
10 Gyr), Vzs are rigidly shifted toward larger values and 
show a steeper trend with stellar mass, with Vz ~ for 
the most massive galaxies. Similarly, Vt. are shifted to- 
ward more negative values while their trend with mass is 
shallower (see T-fTl). Note that in case of no age gradients, 
(negative) shallower Vt.s naturally correspond to (nega- 
tive) shallower colour and metallicity gradients. On the con- 
trary, in presence of non-zero age gradients negative colour 
gradients could correspond to positive M/L gradients due 
to recent star formation episodes. Almost all of these corre- 
lations are significant at more than 99% (see Table [1]). 

Satellite galaxies behave differently from centrals. 
Colour gradients are less dependent on mass, while the V T , 
trend is steeper with mass. This is driven by the V age which 



is increasing at the higher mass bins, being Vz less depen- 
dent on M*. The fact that the Vt, trend is mainly driven 
by the age gradients is confirmed when forcing V agc = 0: in 
this case the trend of Vt, with mass gets shallower although 
Vz show a steep increasing trend with M*. 

In Fig. [2] colour, metallicity, age and M/L gradients 
are plotted as a function of N ga i. We find that in denser 
environments (higher N ga i) central galaxies have shallower 
colour gradients, going from V g _ r ~ —0.09 in the field to 
V g _ r ~ —0.04 in the clusters (N ga[ >50). This trend is 
reproduced by a similar behaviour of metallicity gradients 
which are on average ~ —0.5 in the field and ~ —0.3 in 
denser environments. On the contrary, V age are only slightly 
steeper in densest environments. According to the trends of 
colour and Z gradients, Vt, tend to be steeper in the field 
(~ —0.05), almost null in poor groups (N ga i ~ 10), while 
they are positive (~ 0.1) at N ga i !>50. The trends with en- 
vironment are clearer if we divide the sample in terms of 
the central age, agei. Three age intervals are adopted, i.e. 
< agei ^ 6, 6 < agei ^ 10 and agei > 10, which gather 
1808, 594, 1123 centrals and 577, 198, 366 satellites, respec- 
tively. Older galaxies have progressively shallower age and 
metallicity gradients, independently from N ga i (see also the 
discussion in T+10). The trends with N ga i are still statisti- 
cally significant for Vz, which become shallower going from 
the field to denser environments by ~ 0.15,0.2,0.4, while 
the correlations for V agc get flat for each age bin. The be- 
haviour of Vt, is different: oldest galaxies have negative 
Vt, (~ —0.2, —0.1) with a shallower trend with N ga i (in the 
very rich systems they are almost zero), while the youngest 
ones have positive gradients (~ 0.05, 0.15), which get larger 
in higher density environment. 

On the contrary, satellite galaxies have gradients which 
are completely independent of the environment they live (see 
slopes in Table [2]). To have a complete view of behaviours in 
satellites, we have also considered the rank-2 satellites, i.e. 
the second most massive galaxy in each group and plotted 
their gradients in the right panel of Fig. [2] Here, we see little 
changes in the V g - r and Vz trends, while V age and Vt, 
seem to be steeper functions of N ga i, with steeper values in 
very high density environment. To avoid a too low statistics, 
we have not investigated higher rank systems. 
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Figure 2. Gradients in terms of the environment for centrals (left) and satellites (right). Medians and 25- 75th quantiles are shown. From 
the top to the bottom, colour, metallicity, age and M/L gradients are plotted as a function of N ga i. The black symbols are as in Fig. [l] 
Moreover, blue, green and red symbols are for galaxies with central age in the intervals < age\ ^ 6, 6 < age\ ^ 10 and age\ > 10 Gyr, 
respectively. In the bottom panel, the cyan line is relative to the V-band M/L gradients. Orange lines and yellow regions set medians 
and 25- 75th quantiles for the second most massive galaxy in each group. 



As a comparison, we also show in Fig.[2]the results when 
galaxy colours are fitted to a synthetic spectral model with 
no age gradient. In this case, metallicity gradients turn out 
to be shallower, independently of N ga i, while M/L gradients 
are steeper. Note that these last trends are consistent with 
what was found for the older galaxies in the reference fit. In 
centrals, the trend is almost unchanged for Vz, while a shal- 
lower slope il found for V T » • For satellites, slopes resulted 



unchanged. For completeness, in the bottom panel we also 
compare with V-band M/L gradients, which despite a slight 
shift are very similar to the reference B-band ones. 

We have finally checked that the different trends in Fig. 
[2] between centrals and satellites survive if we cut the most 
massive galaxies with log M*/M© > 11.4 or if we include in 
our sample the galaxies with 10 < log M*/Mq ^ 10.5. 

If we bin the galaxies in terms of stellar mass we find 
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Table 1. Sign of the correlation between gradients and M* or 
N ga i and confidence limits. The significance of the correlation is 
obtained by applying the Student's t-distribution to the Spear- 
man rank factor. 





Centrals 


Satellites 


V E _, - M* 

Vg-r " N gal 


f (99%) 
t (99%) 


T (95%) 
t (95%) 




All free 


Vagc = 


All free 


Vagc = 


Vagc - M. 

V z - M. 
V T>1 - M» 


t (99%) 
t (95%) 
t (99%) 


t (99%) 
t (99%) 


t (99%) 
Null (99%) 
t (99%) 


t (99%) 
t (99%) 


Vagc " N gal 

Vz - N gal 

VX» - Ngal 


t (90%) 
t (99%) 
t (99%) 


t (99%) 
t (99%) 


Null (99%) 
t (95%) 
t (95%) 


t (99%) 
t (95%) 


Table 2. Slopes of the correlation between gradients and M, or 
Ngal an d lo" uncertainties derived by bootstrap method. 




Centrals 


Satellites 


Vg_r - M» 
Vg-r " N ga l 


0.06 ± 0.01 
0.03 ± 0.01 


0.04 ± 0.02 
0.01 ± 0.01 




All free 


Vagc = 


All free 


Vagc = 


Vagc - M, 
V z - M, 
Vt» - M* 


0.21 ± 0.04 
0.23 ± 0.04 
0.20 ± 0.02 


0.30 ± 0.04 
0.12 ± 0.01 


0.28 ± 0.07 
0.17 ± 0.09 
0.26 ± 0.06 


0.26 ± 0.05 
0.10 ± 0.02 


Vagc - Ngal 

V Z - Ngal 
Vl, " Ng^ 


0.08 ± 0.04 
0.12 ± 0.02 
0.11 ± 0.02 


0.13 ± 0.03 
0.05 ± 0.01 


0.02 ± 0.04 
0.05 ± 0.04 
± 0.02 


0.05 ± 0.02 
0.01 ± 0.01 



that if considering the lower mass galaxies the trends in Fig. 
[5] for centrals disappear, while for more massive systems the 
trend become more significant and also the N ga i — M* rela- 
tion is steeper. On the contrary, if we bin the gradients- M* 
relations in Fig.[T]in terms of N ga i we find that field galaxies 
show some trends in V age and Vz, which are stronger in 
group centrals. 



4 DISCUSSION AND CONCLUSIONS 

In this paper we have investigated the correlation between 
the colour, M/L and stellar population gradients with en- 
vironment for a sample of local centr al and satellite SDSS 
central galaxies (|Blanton et al.ll2005bl , . We have found some 
indications which point to a different behaviour in terms of 
stellar mass and the environment for central and satellite 
galaxies. Central galaxies dominate the galaxy sample in 
the massive side, thus have gradients which behave like in 
T+10 and T+ll when the gradients are plotted in terms 
of M*, with Vg-r, Vt» and Vz getting shallower and V agc 
(positive) steeper in the most massive galaxies. On the con- 
trary, satellite galaxies have a V g _ r which is less dependent 
on M*. Thus, while V age are steeper at larger masses (where 
they drive the trend of Vt,), metallicity gradients show a 
shallower trend. The steeper and stronger trend of Vz and 
the slightly shallower Vz in the most massive centrals, when 
compared with satellites, a re consistent with the role of dry 
mergmgs i|Kobavashill2004 ). which are more efficient to mix 
stellar populations at the largest masse s where their frac- 
tion is larger (i.e. Ide Lucia et al. 1120061 ,. Shallower trends 



are found in satellites, suggesting a minor role of such kind 
of phenomena. In satellite we have found a strong trend of 
Va g e with mass (stronger that the one found in centrals), 
with very massive satellites having steeper V ag e, and thus 
younger cores. These surviving young cores could be due to 
bursts of star formation after gas-rich mergers or close en- 
counters. In T+10 we have also seen that central age has an 
important role, since centrally younger galaxies have steeper 
gradients, while older systems have, on average, null V ag e 
and shallower Vz- This seems consistent with the interven- 
tion of processes like dry-mergings, which are responsible for 
a star formation suppression and a mixing of stellar popu- 
latiorQ. 

Plotting the gradients in terms of the number of galax- 
ies inside the group, N ga i, for central galaxies we find that 
age gradients seem less dependent on the environment, 
while colour, M/L and metallicity gradients are shallower in 
denser environment (with a number of galaxies in the groups 
of Ngai ~ 50 or more) when compared with isolated galaxies. 
This correlation is even stronger when galaxies are classified 
in terms of their central ages, which have been shown to 
drive the scatter of the relation between metallicity/age gra- 
dients and mass. On the contrary, massive satellite galaxies 
do not present any trend with N ga i (and independently on 
the age bin adopted). Consistently with the results found 
in terms of stellar mass, while dry-mergings seem to be im- 
portant process in massive centrals, no strong indication is 
found for satellites. Finally, we have also checked that rank- 
ing satellites on the basis of their masses, and considering 
the second most massive galaxy in each group, some slight 
positive correlations seem to appear in our trends in Fig. [2] 
These findings suggest that if on one hand conclusions on 
centrals are firm and clear, on the other hand satellites re- 
sults are more complicated to interpret. More detailed anal- 
ysis would investigate the correlation of population gradi- 
ents not only with N ga i, but also with the local and global 
density, the distance from the center of the groups and the 
mass ranking in the group, which all together give a more 
complete representation of the environmental properties. 

The positive correlation found between colour and Z 
gradients with both Ai» and N ga i in central galaxies are not 
completely independent each others, since it is easy to show 
that more massive central galaxies are residing in denser 
environments with larger cluster/group haloes and popu- 
lated by a larger number of galaxie s fe.g. lWhilev et al.! 2008. 
IStott et al.ll200Sl , IStott et alj|2010h . 

The systematics on the gradients have been investigated 
in the Appendix fS] where we have shown that Vx, is very 
little affected (see, e.g., T+ll), while a shift and a larger 
scatter could be induced in the estimates of V ag c and Vz- 
The trend with M* and N ga i could be possibly affected, in 
the sense that the true trends for Vz can be shallower or 
flat and the one for V age steeper (in the worst cases when 
input data are strongly unaccurate). 



4 A comment is in order here. Metallicity gradients are less sen- 
sitive to recent starburst events, while they are more indicative of 
the total integrated star formation history. This because the metal 
content responds more slowly to the star formation as metals are 
not strongly altered by new stars, but instead by the integral of 
all the star formation history. On the other hand age gradients 
are by definition sensitive to (significant) star burst events. 
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From the observational point of view, these are the 
first steps in the systematic study of the correlations be- 
tween stellar population gradients and the environment 
as measured by N ga i where the behavior of the central 
galaxies is distinguishe d by the one of satellites (see e.g. 
lLa Barbera et all l201ll ). Former studies made no selec- 
tion on the galaxy type (centrals or satellites), and found 
that color gradients in massive gal axies are shallower in 
dense r environment l)Ko fc ~Lmll2005l and 



2005) similarly to our findings, although, 



La Barbera ct al 



La Barbera et al 



201 ll have shown that the optical-IR V g _K is almost 
independent of the env ir onme nt. Along the same line, 
iRoche. Bernardi fc Hyde! (|2010h have demonstrated that 
central galaxies have shallower gradients than normal E/SO 
galaxies at the same luminosity. 

For wha t con cerns stellar population gradients, 
ISpolaor et all |2009l) (their right panel in Figure 1) have 
compared central galaxies in both clusters and groups with 
those in the field, but a clear trend is not found since the 
scatter of the data is too large and the sample too small. 
Using a combined samp le of satellite and central galaxies, 
lLa Barbera et abl (|201 lh have found that V age are slightly 
increasing with M* (similarly to our Fig. while Vz de- 
creases, and galaxies in groups present steeper V age and Vz, 
which are in contrast with our results. The latter discrepant 
results may be a consequence of a different sample selection 
or systematics in our stellar population models. Differences 
in the galaxy sample selection would contribute to some of 
the discrepancies discussed above. In fact, they are adopt- 
ing a sample without any selection in centrals and satellites, 
with the latter dominating in number, while in our case the 
number of centrals is larger and if we derive the plots in 
Fig. [2] for the whole sample, then no trend is found and a 
better agreement is found. Then, if we consider the worst 
case analyzed in simulations discussed in the Appendix (see 
e.g. Fig. IA1|I there seems to be some spurious trends intro- 
duced by the use of the optical bands only: if we correct 
our results for this systematics we can also conclude that 
a steeper trend with mass for V age and a shallower trend 
(if any) for Vz would be accommodated by our data, in a 
way more consistent with lLa Barbera et alj i|201ll ) findings. 
If this is the case, it might indicate a major role of gas-rich 
mergings, cold accretion at high redshift, later gas accretion 
or close encounters, in satellites, in order to produce steeper 
gradients in very massive galaxies, rather than dry merging 
as we would conclude if the gradients are shallower like op- 
tical data suggest. Model systematics seem to weakly affect 
the conclusions with N ga i though, which clearly indicates a 
difference of the gradient behaviour between central galaxies 
and satellites. 

Our results give an indication of the role of physical 
processes in the mass accretion of massive central and satel- 
lite galaxies. The effect of merging and, in general, of inter- 
actions with environment crucially depends on the kind of 
process under analysis (minor vs major merging, dry vs wet 
mergings, accretion, stripping, etc. etc.). E.g., the positive 
age gradients found mainly in both (centrally young) cen- 
tral and satellite galaxies supports two kinds of processes: 
a) a dissipative formation picture, whereby wet mergers fuel 



the central region with cold gas (do Lucia et al 
b) cold accretion at high redshift 



20061). 



Dekel fc Birnbo im 2006), 



which generate the observed younger stellar populations 



in the center. But, central galaxies have more chances to 
merge or interact with neighbors in denser environments, 
in fact they are predicted to have larger masses in clus- 
ters wh en compared with similar galaxies in the centers of 
groups (R omeo et al.ll2008l ) . The flatter gradients of centrals 
we observe in denser environments are consistent with the 
intervention of a larger number of dry (major) galaxy merg- 
ings, which are p redicted to mesh stellar popu lations inside 
the galaxies (e.g. lWhitelll980l : lKobavashill2004l ). On the con- 
trary, satellite galaxies do not show any strong evidence in 
favour of this kind of merging events. The observed younger 
cores in the youngest galaxies suggest that mergings was not 
so effective to erase these pre-existing age gradients (possi- 
bly generated after a gas-rich merging or cold accretion) or 
some recent accretion in the center has happened. While, 
dry mergings were really efficient to dilute those gradients 
in older galaxies. 

The results discussed in this paper represent an useful 
reference exercise for further analysis performed on higher- 
z samples, typically involving optical rest frame bandpasses 
(e.g. from the VLT Survey Telescope, VST). Future analysis 
with wider wavelength baselines will help to alleviate the 
problem of parameter d egeneracies in stellar modelling (e.g., 
lLa Barbera et alJlioTll ). 
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APPENDIX A: 
STELLAR FIT 



SYSTEMATICS IN THE 



Here, we check for systematic effects on our gradient (V agc 
and Vz) estimates from the stellar populati on fits. Because 
of the well known age-metallicity degenera cy (Worthc vl 19941 . 
iBruzual fc Chariot I [20ol . iGallazzi etai~l 120051 ). the stellar 
population parameters and the gradients might be biased 
when using the optical colours only, as we do now. Widen- 
ing the wavelength baseline to include near-infrared (NIR) 
colours should ameliorate the degeneracy. Although the syn- 
thetic prescriptions in the near-IR region of the spectra are 
still uncertain, it is very important t o understand their im - 
pact on our analysis. In T+10 and iTortora et all |2011bT ) 



we checked age and metallicity inferences using optical ver- 
sus optical+NIR constraints, and found, on average, little 
systematic difference. Anyway, depending case by case, a 
spurious shift can appear in the estimated gradients and a 
wide scatter when the fit is made only using optical data. 
Here we will carry out a similar analysis to study the trends 
with stellar mass and N ga i, using a suite of Monte Carlo sim- 
ulations. Because of the large amount of details involved, we 
will consider the results from these simulations as a qualita- 
tive guide to understand in what direction the correlations 
found would change. 

We extract 1000 simulated galaxy spectra from our 
BC03 spectral energy distribution libraries with random, 
uncorrelated stellar parameters and gradients. In order to 
have gradients similar to the ones we observe, we only im- 
pose some constraints on the input stellar parameters, pro- 
ducing, on average, V age ~ 0.2 and Vz ~ —0.3. We add 
simulated (equal) measurement errors for each band, as 
randomly extracted steps from the interval (~S,+S), with 
8 — 0.01, 0.03, 0.05. We apply our fitting procedure and then 
compare the output parameter estimates to the input model 
values. We perform the fit (1) using only the optical SDSS 
bands uqriz and (2 ) adding NIR photometry (J, H and Ks; 
I Jarrett et ail 120031) and (3) u ltraviolet photometry (NUV 
and FUV; iMartin et al] 120051 '). Finally, to plot the results 
in terms of stellar mass we have adopted two linear rela- 
tion to convert the M/Ls into stellar mass. As test relations 
we have recovered from our results at Ri and R2 the fol- 
lowing best relations: log M/L — —0.75 + 0.11 log M* and 
log M/L = -1.08 + 0.15 log Af*. Moreover, the best fitted 
relation between stellar mass and N ga i have been derived 
and used to plot the shift in the gradients in term of N ga i . 

We define the quantity V x,out — Vx,/jv, where X — 
{age, Z), IN and OUT are the input and output X values. In 
Fig. |A1] we show the V x,out — Vxjjv as a function of M» 
and N ga i, for both age and metallicity and the two M/L- 
M* relations. We find that, adding the NIR or UV data, 
the gradients, independently from the shift 8 adopted, are 
perfectly recovered with a very little scatter. With optical 
data only, the uncertainties are larger and some systematic 
shifts emerge, mainly for 8 0.03. In this case, our predicted 
V agc and Vz are underestimated and overestimated with 
respect to the input values, respectively. These discrepancies 
are larger for galaxies with larger M, and N ga j, and in the 
worst case can amount to ~ 0.1 with respect to the ones with 
lower M* and N ga i, for 8 = 0.05. Although these findings 
would led to a possible flattening of the trends of Vz in 
terms of mass and N ga i, a weaker trend would survive. If 
one assumes a constant relation between stellar mass and 
Ng a i, as for satellites, the trend with N ga i in Fig. IA1I would 
be canceled and a simple offset in gradient trends would 
arise in Fig. [2] 
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Figure Al. Systcmatics in stellar population fitting in terms of stellar mass (left panels) and galaxy number in the groups (right panels). 
We define the quantity Vx.OfJT — Vx.IJV: where X = (age, Z), IN and OUT are the input and output values. In each big panel, from 
the top to the bottom we show the systematics in age gradients for two M/L-M* relations and Z gradients for the same two M/L-M* 
relations. From the left to the right we show the cases S = 0.01, 0.03 and 0.05. The continue line is the median while dashed ones are 
the 25- 75th quantiles. Green, red and blue are for the fit using optical, optical+NIR and optical+UV, respectively. 



